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1.0  TECHNICAL  BACKGROUND 

Following  the  development  of  radar  In  the  mid-1940's, 
investigators  began  to  notice  a  pattern  of  enhanced  radar 
reflectivity  near  the  melting  level  in  some  precipitating  clouds. 
This  phenomenon  was  labeled  the  "bright  band"  and  was  the  subject 
of  several  earlier  theoretical  and  observational  investigations 
(for  example,  see  Ryde,  1946;  Cunningham,  1947;  Austin  and  Bemis, 
1950;  Marshall  and  Gunn,  1952;  Wexler,  1955;  Lhermitte  and  Atlas, 
1963).  Battan  (1973)  summarizes  these  earlier  works  which 
concluded  that  the  bright  band  occurred  in  conjunction  with  the 
melting  of  snowflakes  or  an  ice  crystal  aggregation  process  (as 
opposed  to  the  accretional  process  which  forms  graupel  and  hail 
in  convective  clouds).  However,  these  inferences  were  drawn 
without  the  benefit  of  in-si^tu  observations  and  therefore  did 
not  explain  the  physics  of  the  process. 

The  advent  of  PMS  two-dimensional  imaging  probes 
( Kno  1 1  enber  g  ,  1976)  facilitated  the  determination  of  in-si_tu 
hydrometer  size  distribution  from  airborne  measurements  (Gordon 
and  Marwitz,  1984).  This  has  led  to  the  recent  evaluation  of  the 
physical  processes  associated  with  bright  bands  using 
observations  from  aircraft  by  Stewart  et^  a  1_ .  (  1984  ).  This  work 

presented  a  comprehensive  picture  of  thermodynamic  and 
microphysical  processes  associated  with  the  bright  bands  near  the 
Sierra  Nevada  mountains  in  California.  This  report  also  presents 
some  i.n  s_itu  airborne  observations  with  concurrent  radar 
observations  to  help  identify  the  physical  processes  associated 
with  the  melting  layer. 

1.1  Physical  Processes  in  the  Melting  Layer. 

Radar  observations  typically  show  an  enhanced  reflectivity 
(in  bright  band)  on  the  order  of  10-15  dBZ  increase  in 
reflectivity,  located  about  300  m  below  the  0  C  isothermal 
layer.  Several  physical  factors  have  been  proposed  which  could 
theoretically  increase  radar  reflectivity  in  the  bright  band  (see 
Battan,  1973  and  Cohen  and  Sweeney,  1983  for  a  review);  however, 
an  evaluation  of  the  most  recent  data  strongly  suggests  that  a 
shift  in  particle  size  distribution  toward  larger  sizes  is  the 
most  probable  explanation.  Since  the  aggregation  efficiency  is 
directly  related  to  temperature  until  melting  occurs,  and  since 
recent  observations  by  Stewart  et  aj_^_  (  1984  )  show  a  strong 
correlation  between  aggregation  and  the  bright  band  effect,  there 
is  strong  evidence  that  the  bright  band  is  associated  with 
aggregation.  Figure  1.1  is  a  reproduction  from  Stewart  et  al^ 
(1984)  showing  the  correlation  between  ice  particle  size 
(measured  by  the  PMS  2-D  probes)  and  radar  reflectivity.  They 
also  point  out  that  another  factor  which  may  support  the 
aggregation  process  is  an  ice  multiplication  process  active  at 
-5°C  (Hallett  and  Mossop,  1974).  As  shown  in  Figure  1.1,  there 
exists  a  threefold  increase  in  ice  crystals  (needle  shaped)  near 
- 5 ° C .  These  needles  may  form  aggregates  or  combine  with 
existing  dendritic  crystals  to  form  large  lacy  aggregates.  These 
large  aggregates  would  presumably  tend  to  clump  together  near  the 
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Simultaneous  profiles  of  reflectivity  factor  z  and  root-meansquare  particle 
fall  speed  in  light  (I  mm/hr).  Steady  precipitation  with  a  bright  band.  From  Lhermitie  and 
Atlas  (19611. 


Z«.  «e(Zt) 

Vertical  profile  of  radar  reflectivity,  and  the  cross  section 
of  the  King  Air  flight  track  from  0711  PST  (takeoff)  to  0725  PST 
on  2  March  197$.  The  radar  data  were  computed  from  a  7.5* 
elevation  scan  of  the  CP-3  radar  at  0719  PST  and  from  a  12-5* 
elevation  angle  scan  at  0919  PST.  (From  Stewart  £t  al.  (1984) 


Figure  1.1  Examples  of  Bright  Band  Radar  Reflectivity 
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melting  layer  and  melt  into  large  drops  which  enhance  radar 
reflectivity. 


The  effect  of  water-coated  ice  particles  on  radar 
reflectivity  has  been  studied,  and.  according  to  Ekpenyong  and 
Srivastava,  (1970),  when  one-third  of  the  original  volume  of  ice 
has  melted  into  a  water  film  surrounding  the  ice,  the  index  of 
refraction  will  have  changed  from  ice  to  that  of  water.  Cohen 
and  Sweeney  estimate  that  this  would  result  in  a  33  dBZ  increase 
in  radar  reflectivity  (much  larger  than  the  observed  effect). 
However,  Knight  (1979)  found  that  a  circular  shell  of  water  only 
forms  in  the  advanced  stage  of  melting.  Consequently,  the 
possibility  exists  that  an  effective  increase  in  hydrometeor 
diameter  resulting  from  a  less  dense  ice  center  could  enhance 
radar  reflectivity.  This  process  is  difficult  to  ascertain  from 
aircraft  observations. 


Another  feature  associated  with  bright  bands  is  the  relative 
decrease  in  radar  reflectivity  below  the  melting  layer. 

The  increase  in  fall  speed  of  the  rain  drops  below  the  melting 
layer  (see  Figure  1.2)  will  create  an  effective  decrease  in 
number  concentration.  Lhermitte  and  Atlas  (1963)  found  that  this 
accounted  for  about  one-third  of  the  observed  decrease  in 
reflectivity  below  the  melting  layer  (about  8  dBZ)  was 
accounted  for  by  the  increase  in  fall  velocity  (4  dBZ)  and  the 
change  in  shape  factor  to  spherical  (about  3  dBZ).  Stewart 
et_^  li  •  also  concluded  that  raindrop  collisional  breakup  below 
the  melting  layer  was  significant  and  tended  to  shift  the  drop 
size  spectra  toward  a  Marshall-Palmer  distribution  with  the 
increasing  temperature.  The  decrease  in  observed  radar 
reflectivity  was  in  agreement  with  the  shift  in  the  size  spectra 
toward  smaller  diameters. 


The  aforementioned  observations  discussed  in  the  literature 
suggest  that  several  physical  processes  which  affect  the  melting 
level  (and  the  ensuing  bright  band)  are  complex  and  interactive. 
It  is  also  apparent  that  the  physical  processes  which  affect  the 
melting  layer  have  their  origins  at  colder  temperatures  ( i . e . , 
ice  multiplication  and  aggregation)  and  extend  well  below  the 
bright  band  (i.e.,  increased  terminal  velocity  and  collision 
breakup).  In  the  present  work  we  discuss  airborne  measurements 
in  the  melting  layer  in  terms  of  this  chain  of  events. 


Overview  of  Airborne  Measurement  Systems 


For  this  study  of  the  melting  layer  CIC  provided  the  PMS 
twin  engine  t  u  r  b  o - c  h  a  r  g  e  d  Beech  Baron  (Figure  1.3)  which  is  fully 
de-iced  and  has  a  good  performance  capability  at  the  altitudes 
covered  in  this  project.  (See  Section  1.3  for  the  Beech  Baron 

aircraft  specifications. ) 


The  Baron  carries ,  in  standard  research  configuration,  a 
full  compliment  of  PMS  partic’e  sizing  and  imaging  probes,  a  gust 
boom  and  the  conventional  state  parameter  sensors  A  multiple 
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Schematic  profile  illustrating  significant  trends  of  various  microphysical  and  state 
parameters  through  a  typical  melting  layer  as  deduced  by  synthesizing  the  observations  of  all 
the  case  studies. 


Figure  1.2  Profile  Illustrating  Trends  of 

Microphysical  and  State  Parameters 
at  the  Melting  Layer  (from  Stewart 
jjt  al .  ,  1984) . 
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DHE  position  keeping  system  was  installed  which  cycles  through 
six  VORTAC  stations  every  six  seconds,  and  provides  position 
accuracy  of  better  than  0.5  km  in  straight  flight.  A  forward 
looking  Ka  band  radar  was  installed  in  the  nose  of  the  Beechcraft 
Baron  (in  place  of  the  standard  radar).  The  return  video  was 
processed  by  fifteen  range  gates,  displayed  in  real  time  on  the 
cockpit  display,  and  recorded  on  magnetic  tape. 

The  data  were  recorded  on  compu t e r - compa t 1 b 1 e  nine  track 
magnetic  tape.  A  CRT  display  presented  real-time  data  in 
engineering  units  to  the  scientist  in  the  cockpit.  A  hard-copy 
printer  output  of  the  data  was  also  available  to  the  scientist  in 
flight. 

The  post  flight  data  processing  was  done  on  the  CIC  IBM-XT 
and  the  AFGL  Vax  systems.  The  AFGL  system  also  benefitted  from 
the  use  of  several  man-years  of  data  analysis  software  developed 
by  CIC  and  the  University  of  Wyoming  during  the  HIPLEX  project 
(see  Appendix  B).  Examples  of  the  computer  products  are  given  in 
Section  3.0. 

1.3  Performance  of  Beech  Baron  Aircraft 

The  performance  characteristics  of  the  Beech  Baron 
research  aircraft  are  shown  in  Table  1.1. 


Table  1 . 1 

TURBOCHARGED  BEECHCRAFT  BARON  MODEL  58-TQ  N23712 


Power,  each  engine 

310  hp 

Maximum  speed 

235  kts. 

Fuel  capacity 

175  gal 

Service  ceiling 

25,000  ft 

Single  engine  service  ceiling 

13,000  ft 

Climb  at  gross  wt 

1500  ft/min 

Maximum  endurance 

5  hrs 

Maximum  range 

900  nm 

Take-off  distance 

2400  ft 

Landing  distance 

2500  ft 

Gross  weight 

6100 

Empty  weight 

4570 

Seating  (Restricted  Catagory) 

3  crew 

Length 

30  ft 

Wing  span 

38  ft 

1.4  Cloud  Physics  Instrumentation. 

1.4.1  Particle  Spectrometers. 

The  PMS  research  aircraft  is  equipped  with  a  Forward 
Scattering  Spectrometer  Probe  (FSSP),  Optical  Array  Two 
Dimensional  Cloud  Particle  Imaging  Probe  (2DC)  and  a  Two 
Dimensional  Precipitation  Imaging  Probe  (21)  P )  .  The  FSSP  sizes 
and  counts  particles  from  0.5  to  45  microns,  the  2  D  -  C  sizes  and 
counts  particles  from  25  200  microns  and  the  21)  P  operates  in  the 
range  from  200-6400  microns.  Figure  1.3  shows  the  PMS  particle 

spectrometer  probes  mounted  on  the  wing  tips  of  the  Baron.  C1C 
has  previously  developed  copious  software  to  process  PMS 
spectrometer  probe  data;  samples  of  computer  products  are  shown 
In  Section  3. 

The  PMS  2D-C  probe  is  equipped  with  a  depolarization  option 
which  helps  to  discriminate  ice  and  water  particles.  The 
depolarization  feature  works  on  the  principle  that  a  birefringent 
substance  such  as  ice  will  change  the  polarization  vector  of 
incident  polarized  light,  whereas  water  drops  will  not 
significantly  affect  the  polarization  vector.  Although  the 
discrimination  is  not  absolute  in  all  cases  ( i . e . ,  large 
asymmetrical  drops  can  cause  a  change  in  the  polarization  vector) 
the  signal  is  a  good  indication  of  the  relative  phase  of 
hydrometeors  (Lawson  and  Stewart.  1983). 

1.4.2  State  Parameter  and  Air  Motion  Sensing 
Instrumentation. 

Table  1.2  lists  the  performance  characteristics  of  the 
research  instrumentation  used  by  C I C .  All  sensors  and  systems 
were  calibrated  at  the  Cl C /PMS  facility  prior  to  the  field 
season.  A  tower  flyby  was  conducted  at  the  instrumented 
BAO/NOAA  tower  located  10  miles  east  of  the  Boulder  Airport  to 
verify  the  dynamic  calibration  of  pressure,  temperature  and  dew 
point  temperature  prior  to  the  field  season. 

1.5  Ka  Band  Weather  Radar  , 

CIC  provided  a  Westinghouse  AN/ A  PI)  7  Ka  band  radar  for  use  on 
this  project  (see  Table  1.3  for  radar  characteristics).  The 
existing  X-Band  radar  was  removed  and  the  API)  7  radar  was  located 
in  its  place. 

Radar  signal  processing  was  accomplished  using  a  n 
integrating  sample  and  hold  circuit  (  r  e  f  e  r  r  e  <1  to  as  a  v  i  d  e  o 
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TABLE  1 . 3 


WESTINCHOUSE  AN/APD-7  Ka  BAND  WEATHER  RADAR  SPECIFICATIONS 


R-F  Operating  Frequency - 

Peak  Power  Output - 

Average  Power  Output - 

Pulse  Repetition  Frequency 

Pulse  Width - 

Antenna  size - 

Beat*  width - 

Minimum  detectable  signal- 

Range  gate  size - 


34 . 9  GHz 


100  KW 


14  watts 

1845  pr f 

0.1  usee 

12"  dia. 

1.5° 

8  dBZ  at  1  km 


15  gates :  600 ,  1800 
3000.4200,5400m;  6.6 
7.8,  10.2,  12.6,  15, 


integrating  processor  or  VIP)  which  had  previously  been  used  by 
CIC  on  several  of  its  ground  based  radar  systems. 

The  radar  transmitter  generated  a  0.1  micro  second  pulse  and  the 
return  video  was  integrated  and  recorded  every  1.0  seconds.  A 
series  of  range  gates  were  recorded  with  the  first  gate  set  at  600m 
ahead  of  the  aircraft,  and  subsequent  gates  spaced  as  shown  in 
Table  1.3. 

16  Navigation  Systems. 

1.6.1  Multi-DME  Position  System. 

CIC  developed  a  multi-DME  position  system  for  use  in 
positioning  the  research  aircraft.  The  system  recorded  distance 
data  from  up  to  six  different  ground-based  VORTAC  stations  at  a 
rate  of  one  station  per  second.  The  KDM-7000  DME  used  in  this 
system  has  a  resolution  of  20  m.  and  the  system  accuracy  is  on 
the  order  of  200  m.  when  four  or  more  stations  are  received 
(Lawson  et^  a  J.  .  ,  1980).  An  additional  advantage  was  that  there 

is  no  drift  with  time  (as  in  INS  positioning  systems)  since  the 
absolute  data  are  refreshed  every  six  seconds.  The  area  around 
Hanscom  AFB  is  densely  populated  with  VORTAC  stations  as  shown  in 
Figure  1.4.  At  normal  flight  levels  around  Hanscom  AFB  DME 
signal  reception  was  adaquate  and  gave  good  positions  for  this 
melting  layer  project. 

1.6.2  VOR/DME. 

The  CIC  data  system  routinely  recorded  aircraft  VOR  and 
DME  data  on  the  research  data  system,  as  well  as  multi-DME. 

Position  from  VOR/DME  is  only  accurate  to  about  1  km. ,  but  the 
data  are  very  reliable  and  provided  useful  redundancy  (the  multi 
DME  system  failed  on  some  days)  . 

1.7  Post  Season  Data  Processing. 

CIC  copied  all  the  raw  data  tapes  onto  9 -track  1600  bpi 
computer  compatible  magnetic  tape  and  delivered  these  tapes, 
along  with  all  data  reduction  programs  (written  in  Fortran  and 
formatted  in  ASC  II  on  9-track  computer  compatible  magnetic  tape) 
to  AFGL  personnel.  All  notes  recorded  by  the  on-board 
scientist  were  delivered  to  AFGL.  CIC  and  OPHIR  personnel  then 
continued  to  work  with  AFGL  personnel  to  develop  analysis 
software  and  techniques  to  elucidate  the  microphysical  and 
thermodynamic  processes  operating  within  and  near  the  melting 
layer  . 

1  . 8  AFGL  M  -meter . 

The  AFGL  M  meter  estimates  the  integrated  mass  of 
hydrometeors  which  are  impacted  on  a  rapidly  spinning  disk  by 
measuring  the  change  in  angular  momentum  of  the  disk.  The  M 
meter  was  installed  on  the  starboard  wing  pylon  which  also 
carries  the  F’MS  FSSP  probe  (see  Figure  1  3).  Proper  wiring  was 


1  1 


radius  of  Hanscom 


already  available  to  provide  electrical  power  and  to  transfer  the 
output  signal  directly  to  the  PMS  Data  Acquisition  System  (DAS). 


In-flight  data  were  available  on  the  cockpit  CRT  and  the 
real-time  printer.  Thus,  fie’d  evaluation  of  the  unit  was 
conducted  on  a  flight-by-flight  basis.  The  primary  evaluation 
tool  was  comparison  of  the  unit  output  with  the  PMS  probes, 
however,  lack  of  sufficient  funding  preempted  any  detailed 
analysis,  and  inquiries  on  this  matter  should  be  directed  to 
Mr.  Vern  Plank  of  the  AFGL. 


2.0  FIELD  OPERATIONS. 

2.1  Period  of  Operations. 

The  aircraft  operated  from  Hanscom  Field  from  March  17,  1986 

through  May  22,  1986.  The  aircraft  was  stationed  at  the  Piper 

fixed  base  operator  at  Hanscom  Field. 

2.2  Air  Traffic  Control  Procedures. 

Prior  to  the  field  program  CIC's  pilot  visited  the  Boston 
ATC  center  in  order  to  explain  the  proposed  flight  operations  and 
to  obtain  their  cooperation  in  filing  flight  plans.  The  aircraft 
was  given  a  code  name  "Cloudy  I."  This  code  name  helped  to 
identify  to  the  FAA  controllers  the  special  aspects  of  this 
flight  opera  t i on . 

2.3  Flight  Operations. 

Table  2.1  is  the  flight  summary  for  this  project.  A  total 
of  53.5  hours  were  flown  on  test  flights,  ferry  to  and  from 
Hanscom  and  on  observational  flights. 

A  standard  flight  profile  was  established  for  the  purpose  of 
measuring  the  chain  of  events  of  the  cloud  physical  passes  as 
they  occur  above,  through,  and  below  the  0  C  level. 

Figure  2.1  shows  an  example  of  the  profile  chosen.  After  take¬ 
off  the  aircraft  ascended  to  the  -4  C  level.  Then,  in  a 
series  of  passes,  the  aircraft  would  alter  its  altitude  to 
pass  through  the  vertical  depth  of  the  melting  layer. 

The  horizontal  tracks  were  designed  such  that  the  center  of 
the  vertical  profile  would  occur  over  Sudbury  where  ground-based 
precipitation  sensors  and  weather  radars  were  stationed 
However,  due  to  air  traffic  control  problems,  not  all  profiles 
could  be  accomplished  as  planned. 


2.4  Data  Tapes. 

All  data  tapes  were  copied  and  given  to  the  Air  Force 
personnel  within  48  hours  after  eacli  flight  The  observer  flight 
notes  were  also  provided.  Copies  of  these  notes  were  also 
retained  by  CIC  for  further  post  field  analysis 
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TABLE  2  .  1 

CIC  FLIGHT  SUMMARY 


PI# 

Date 

Clock  hrs 

Pros  Total 

Type 

PRE  1 

3-15 

1  .  2 

1  .  2 

CA1-IU5 

PRE  2 

3-19 

0 . 5 

1  .  7 

CAL-KC 

FERRY 

3-16-20 

10.3 

13.0 

I V5-BED 

01 

3-22 

0 . 7 

13.7 

FAM-BED 

02* 

3-27 

0 . 8 

14.5 

TEST-ABORT 

03* 

3-27 

1  .  6 

16.1 

TEST-DATA 

04* 

4-7 

2 . 5 

18.6 

DATA 

05 

4-8 

2 . 3 

20 . 9 

DATA 

06* 

4-17 

3 . 4 

24 . 3 

DATA 

07 

4-21 

1  .  5 

25 . 8 

DATA/ABORT 

08* 

4-26 

3 . 6 

29 . 4 

DATA 

09* 

4-27 

3 . 6 

33 . 0 

DATA 

10 

5-2 

0 . 6 

33 . 6 

M-M  TEST 

11* 

5-12 

2 . 3 

35 . 9 

DATA 

12 

5-14 

0 . 9 

36 . 8 

M-M  TEST 

13 

5-14 

1  .  5 

38.3 

M-M  TEST 

14* 

5-17 

2  .  2 

40.5 

DATA 

15* 

5-22 

2  .  1 

42.6 

DATA 

16 

5-22 

1  .  7 

47.3 

FERRY 

17 

5-23 

3 . 5 

50.8 

FERRY/DATA 

1  8 

5-23 

2  .  7 

53.5 

FERRY 

*  Flights  which  provided  data  of  particular  interest 
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3.0  SOFTWARE  FOR  DATA  REDUCTION  AND  ANALYSIS. 


3 . 1  Background . 

The  instrumented  aircraft  collected  data  on  magnetic  tape 
which  was  transcribed  to  hardcopy  computer  products  after  each 
flight.  These  computer  products  provided  the  AFGI,  with  a  quick - 
turnaround  evaluation  of  the  mission  performance  (for  example  see 
Figure  3.1). 

During  a  meeting  on  13  May  1986  with  AFGL  and  CIC  personnel, 
it  was  agreed  that  CIC  personnel  would  pursue  the  analysis  of 
data  after  the  close  of  the  field  season  (an  extension  of  the 
field  season  by  the  AFGL  due  to  lack  of  suitable  weather  reduced 
the  original  estimate  of  funds  available  for  data  analysis). 

AGFL  personnel  recommended  that  CIC  develop  software  to  produce  the 
following  computer  products: 

1.  Plan  view  of  aircraft  position  as  a  function  of  time. 

2.  Integrated  total  water  content  and  equivalent  radar 
reflectivity  including  artifact  rejection  from  the 
PMS  FSSP,  2DC  and  2DP  particle  imaging  data. 

3.  Time  series  plots  of  the  airborne  Ka-band  radar 
reflectivity  data  which  have  been  t i me - c o r r e 1  a t ed 
for  distance  from  the  aircraft. 

4.  Spectral  analysis  of  correlation  between  the  airborne 
Ka-band  radar  data  with  particle  image  data. 

It  was  further  agreed  that  CIC  would  use  these  software 
tools  to  produce  plan-view  and  vertical  cross-section  analyses  of 
the  melting  layer  data  identifying  the  location  and  type  of 
hydrometeors,  the  depth  of  the  melting  layer  and  the  precipitation 
processes  which  were  active.  These  j_n^si_tu  data  were  intended  to 
be  correlated  with  remotesensing  data  from  ground-based  radars 
and  satellite  beacons. 

3.2  Computer  Software  Products. 

3.2.1  PMS  2D  Images. 

Figure  3.2  is  an  example  of  PMS  2DC  and  2 DP  image  data 
collected  on  27  March  1986.  The  images  are  underscored  by  lines 
of  meteorological  data  collected  by  the  aircraft  at  the  same 
time. 

3.2.2  Aircraft  Track 

Figure  3.3  is  an  example  of  a  plot  of  aircraft  t  r a  r  k 
from  the  multiple  D  M  E  system  as  a  function  of  time.  Since  the 
aircraft  flew  repeated  legs  over  the  same  location  on  the  ground, 
it  is  necessary  to  separate  a  mission  into  several  time  periods 
for  the  sake  of  clarity.  Appendix  A  contains  the  flight  tracks 


AIR  FORCE  GEOPHYSICS  LABORATORY  MELTING  LAYER  DATA  COLLECTED  ON  3/27/86  BY  CIC  CLOUD  PHYSICS  AIRCRAFT  N23712 


TIME 

PRES 

ALT 

VOR/OME 

84-DHE 

STA 

PITCH  HOG 

IAS 

TAS 

RSA4T 

RFT 

1034  0 

929 

721 

313/ 

0 
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Figure  3.2  Example  of  images  from  I'.VS  ?I’-C  |>r<>be  (  !  at  i-t-r  images  with 

vertical  scale  of  800  mi  crons)  and  211-1’  prebe  (smaller  images 
with  vertical  scale  of  *> .  mm). 


1 ') 


MLn.  AA  * 


ft. ' 


•  .  -  *  '  >  -  •  '  •  •  -  -  • 

V <*.N  A  fcV^*V  »W 


*r.' 

,  M  «l» 


/V' 

■is 


:>S 


y./i 

V. J 
V. . 
V/J 


I 


-  a  i 


.  .N  ,N  "k  ,%  .S  >  %  ■-  .m 


for  the  priority  days. 

3.2.3  PMS  Hydrometeor  &  Ka-band  Radar  Data  Comparison. 

Figures  3.4  -  3.5  are  examples  of  the  correlation 
between  the  reflectivities  from  the  airborne  Ka-band  radar  and 
the  equivalent  radar  reflectivity  integrated  from  the  PMS 
hydrometeor  spectra.  The  hydro meteor  data  are  corrected  for 
artifact  rejection  (artifacts  created  by  particles  colliding  with 
the  probe  tips,  cloud  water  streaking  across  the  mirrors,  etc.) 
using  the  algorithm  reported  by  Cooper  (1978).  The  airborne  Ka- 
band  radar  sees  hydrometeors  well  ahead  of  the  aircraft  (from  0.6 
to  36  km),  and  these  data  are  adjusted  in  time  using  the  recorded 
true  airspeed  to  correspond  with  the  PMS  particle  imaging  data. 
The  data  in  Figures  3.4  &  3.5  show  a  casual  correlation  between 
equivalent  reflectivity  from  the  integrated  hydrometeor  spectra 
and  the  first  range  bin  of  the  airborne  data,  but  we  fail  to  see 
any  correlation  extending  beyond  the  first  range  gate.  We  feel 
that  this  is  due  to  the  very  light  precipitation  (maximum  of  2  mm 
diameter  drops),  and  that  the  reflected  radar  signal  was  lost  in 
the  noise  past  the  first  range  gate. 

Figure  3.6  shows  an  ensemble  plot  of  equivalent  radar 
reflectivity  from  the  2DP  probe  versus  the  first  range  gate  of 
the  airborne  radar,  and  indicates  that  the  data  contain  a  large 
degree  of  variance. 

3.2.4  Spectral  Analysis. 

The  correlation  between  two  variables  can  be  examined 
in  the  frequency  domain  using  FFT  analysis  techniques.  Figures 
3.7a,  b  and  c  show  the  correlation  between  the  first  range  gate 
from  the  airborne  Ka-band  radar  reflectivity  (adjusted  in  time) 
and  the  equivalent  reflectivity  from  the  integrated  spectra  seen 
by  the  2DP  probe.  The  coherence  plot  shows  significant 
correlation  out  to  about  0.05  Hz,  which  corresponds  to  1.6  km  at 
a  true  airspeed  of  80  m/s.  This  is  some  indication  of  the  degree 
of  stationarity  in  the  precipitation  field  during  this  mission. 
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Figure  3.4  Example  of  time  series  data  showing  equivalent  reflectivity 
from  PMS  2D-P  image  data  (top)  and  reflectivities  from  the 
first  three  range  gates  of  the  airborne  Ka-band  radar  which 
have  been  shifted  in  time  to  correspond  with  the  position 
in  space  of  the  image  data  collected  on  5/12/86. 
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40  CASE  STUDIES 

4.1  Overview. 

As  pointed  out  in  Section  10.  understanding  the  physical 
chain  of  events  leading  to  the  formation  of  precipitation  is  the 
major  concern  in  describing  characteristics  of  the  melting  layer 
During  the  field  observational  period  the  Karon  research  aircraft 
collected  data  at  various  levels  above,  below  and  near  the  0  C 
level.  The  purpose  of  this  section  is  to  examine  the  images  as 
seen  by  the  PMS  2D  probes  and  correlate  these  observations  with 
the  formation  of  precipitation. 

Vertical  cross  sections  of  the  2  DP  data  have  been  analyzed 
and  for  each  leg  of  the  flight  track  the  2  DP  images  have  been 
studied  to  determine  the  occurrence  of  various  types  of  solid 
precipitation  and  to  estimate  the  elevation  at  which  these 
hydrometeors  melted  Since  funding  for  analyses  of  these  cases 
is  limited,  we  have  chosen  to  present  a  brief  description  of  the 
five  cases  selected,  with  a  more  detailed  analysis  of  a  portion 
of  the  27  April  c  ase  T fi  i  s  case  provided  the  best  example  of  a 
steady  state  episode1  which  allowed  us  to  draw  some 
generalizations  The  other  cases  generally  support  the  major 
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especially  near  the  western  end  of  the  track. 

The  aircraft  next  climbed  to  2  16  0  m  (2.2°C)  and  started  leg 
four  in  a  NE  direction.  On  this  leg  there  were  large  u n melted 
aggregates  mixed  with  small  raindrops.  Some  of  the  aggregates 
were  quite  large.  The  aircraft  then  ascended  to  about  2200  m 
(1.4  C)  and  started  leg  five.  Again  on  this  leg  there  were  large 
unmelted  aggregates  mixed  with  small  raindroos. 

Leg  six  was  flown  at  2300  m  where  the  temperature  was  just 
about  0  C.  The  images  were  a  mixture  of  medium  sized  aggregates 
and  graupel  particles,  which  did  not  co-exist  but  appeared  to 
occur  in  patches,  a  re-occurring  theme  found  in  the  data.  There 
was  a  large  variation  in  particle  concentration  along  this  leg. 

Leg  seven  was  flown  at  2400  m  (  0.4  C  and  had  images  very  similar 
to  leg  six  (mixture  of  aggregates  and  graupel).  Leg  eight  was  a 
reciprocal  leg  of  seven  and  showed  the  same  types  of  images  as 
did  leg  nine. 

Leg  ten  started  at  2500  m.  then  the  aircraft  ascended  to 
about  2700  m  (-1.2  C).  The  images  were  a  mixture  of  large 
aggregates  and  small  graupel  The  2DC  malfunctioned  and  the  image 
data  were  unusable  on  this  leg.  Next  the  aircraft  ascended  to 
2800  m  (-1.8  C)  and  made  legs  eleven  and  twelve  at  this 
altitude.  The  images  on  these  legs  were  mostly  graupel  partlc'es 
with  a  few  large  aggregates,  and  some  needles 

At  0342  EST  the  aircraft  descended  to  about  1800  m  (2.9°C  and 
conducted  leg  13.  On  this  leg  there  were  mostly  1  3  mm  raindrops 

with  a  few  Large,  unmelted  aggregates. 

4.4  Case  Study  III  (427  86). 

This  case  has  the  most  extensive  set  of  cloud  physics  data 
with  25  legs  flown  in  three  separate  vertical  profiles.  Table 
4.1  describes  the  first  11  legs  on  this  flight.  The  aircraft 
ascended  from  1800  to  2700  m  while  making  transects  from  Hans  com 
A F B  to  Sudbury  The  principle  precipitation  mechanism  during 
this  profile  seemed  to  be  an  aggregation  process.  Large 
aggregates  (2  5  mm)  were  observed  to  form  just  above  the  0  C 
level  and  melted  forming  1  2  mm  rain  after  descending  about  800  m 

Occasional  large;  melted  aggregates  were  observed  along  w  i  t  h 
t.  he  rain  drops.  There  were  few  graupel  particles  observed  during 
this  profile,  and  the  accretional  growth  of  ice  into  graupel 
part i r I ps  did  not  appear  to  be  a  significant  precipitation 
mechanism  This  vertical  profile  is  notably  different  from  th.it 
observed  later  (0120  0214)  w h i c  h  showed  discrete  pockets  of 

either  graupel  or  aggregates 
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TABLE  4.1.  Description  of  Hydrometeors  Observed  by  the  PMS 
Probes  on  the  Baron  Aircraft  During  the  First 
Vertical  Profile  on  27  April  1986 


L  e  ^ 

IloieiEsii 

A  1 1  i  t  u  <1  e  (  m  ) 

Image  Characteristics 

1 

0014-0019 

1 800  (  2  .  0°C ) 

mostly  rain,  few  large 
unmelted  aggregates 

2 

0019-0024 

1  8  0  0  ( 2 . 0  °  C  ) 

same  as  above 

3 

0025-0028 

1 800  ( 2 . 0°C ) 

same  as  above 

4 

0028-0032 

2135  (0. 3°C) 

mostly  rain,  frequent 
large  aggregates 

5 

0033-0038 

2135(0. 3°C  ) 

same  as  above 

6 

0038-0048 

2275 ( 0 . 4  °C ) 

pockets  of  large  aggregates, 
small  stellars ,  small  drops 

7 

0048-0053 

2275(0. 4  °  C ) 

aggregates,  ice  particles 

8 

0055-0100 

2431(  0. 6  °  C ) 

large  aggregates,  few  small 
rimed  ice  or  graupel  embryos 

9 

0101  -0106 

243  1  (  -0  .  6°C ) 

mixture  aggregates, 
graupel  embryos,  stellar  ice 

1  0 

0108-0113 

2  6  0  0  ( -  1  . 3  °  (  ) 

many  small  aggregates  made 
of  needles,  ice  particles 

1  1 

01  14-01 19 

2  5  8  6  (  1  .  2  °  (  ) 

few  large  aggregates, 

mostly  small  Ice,  some  graupel 

embryos . 


«*, 

A 


« 

l 


» mW.»  »v»  * »  * 


Table  4.2  lists  the  characteristics  of  the  next  eleven 
legs.  This  profile  characterizes  the  pattern  which  we  haze  seen 
in  several  other  portions  of  the  data,  in  that  graupel  particles 
and  aggregates  appear  to  exist  in  separate  pockets  across  the 
cross  section.  The  mechanisms  responsible  for  the  generation  of 
these  inhomogeneous  regions  are  beyond  the  scope  of  this  analysis 
(the  aircraft  did  not  ascend  above  -  3  C  .  where  the  formation  of 
the  graupel  and  ice  particles  originated). 

Figure  4.1  is  an  analysis  of  the  vertical  cross  section  of 
the  data  from  0120  0214.  which  suggests  that  there  were  two 

distinct  regions  of  hydrometeor  types.  In  the  southwestern  region 
(near  Sudbury),  relatively  high  concentrations  (up  to  30  1,  of 
ice  and  graupel  particles  were  observed  above  the  melting  layer, 
while  over  the  northeastern  area  (near  llanscom  AFB),  smaller 
concentrations  of  ice  particles  and  aggregates  were  found  above 
the  0  °  C  level.  Examples  of  the  images  from  these  two  regions 
are  shown  in  Figure  4.2.  The  graupel  particles  co  -  existed  with  a 
significant  region  of  supercooled  liquid  water,  while  the 
aggregates  did  not.  This  regime  seemed  to  persist  in  a  spatial 
relationship  commensurate  with  the  shear  of  the  mean  wind  as  the 
aircraft  ascended,  and  therefore  appears  to  be  a  quasi  steady 
state  feature.  This  is  substantiated  by  the  rapid  spiral  descent 
from  2700  -  1800  m  at  2115  which  showed  the  evolution  of  graupel 

to  small  rain.  Below  the  melting  layer,  the  2DP  images  suggest 
that  the  graupel  melted  into  fairly  uniform,  small  (1  mm) 
raindrops  at  a  level  very  close  to  0  C.  On  the  other  hand,  the 
aggregates  did  not  melt  uniformly  and  some  partially  melted 
aggregates  are  seen  down  to  +3  C ,  where  there  is  a  larger 
variance  in  the  precipitation  spectrum  and  several  larger  (1-3 
mm)  raindrops  are  observed.  Examples  of  images  from  these  two 
regions  below  the  melting  layer  are  also  shown  in  Figure  4.2.  It 
thus  appears  that  the  aggregation  process  may  have  been  more 
efficient  In  utilizing  the  available  supercooled  cloud  water, 
generating  larger  particles  which  took  longer  to  melt  and 
subsequently  produced  larger  raindrops.  This  is  contrary  to 
observations  in  summertime  cumuli  (Cooper  and  Lawson,  1984)  and 
is  discussed  further  in  Section  4  .  6  . 
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TABLE  4.2.  Description  of  Hydrometeors  Observed  by  the  PMS 
Probes  on  the  Baron  Aircraft  During  the  Second 
Vertical  Profile  on  27  April  1986. 


Le£  liieiESTl 


1  0120-0124 


0124-0129 


A  l_t  i_t  udc 


1800(3.7C) 


1800(3.7  C) 


^)}3racteri_sti_(:s 

1  mm  rain  SW; 

1-2  mm  rain  &  unmelted 
aggregate  NE. 

same  as  above 


0130  -0135 


2130(1.8  C ) 


same  as  above 


0136-0141 


2130(1 ,8C) 


same  as  above 


0141-0146 


6  0146-0152 

7  0153-0158 


8  0158-0204 


9  0204  0209 


10  0209-0214 


2244(0.7  C) 


2275  (0.3  C) 
2  4  4  5  ( 0  °  C  ) 


2400  (-0.6  C) 


2650  (  -1  .  5  C  ) 


2757  (  2 . 0  C ) 


11  0215-0220(Spiral  Descent  over  SDD) 


mixed  1  mm  graupel 
&  rain  SW ; 
mixed  1-2  mm  rain 
&  2-5  mm  aggregate  NE 

same  as  above 

same  as  above 

ice  &  0.5  1  mm  graupel  SW; 

ice  and  25  mm  aggregate  NE 

ice  and  0.5  mm  graupel  SW, 

ice  and  1  mm  aggregate  NE 

same  as  above 

Graupel  increasing  from 
0.5  to  2  mm  and  melting 
to  form  1  mm  rain 
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4.5  Case  Study  IV  (5-17-86) 


The  flight  profile  for  this  case  was  very  similar  to  the 


previous  cases.  There  were  six  significant  periods  during  this 
flight  and  they  are  summarized  and  shown  in  Table  4.3. 

TABLE  4.3.  Description  of  Hydrometeors  Observed  by  the  PMS 
Probes  on  the  Baron  Aircraft  on  17  May  1986. 

Leg  lilS§iIST2  Ternge ra_t u r e Xm Image  Characteristics 

1  0300-0310  0.4°C  large  unmelted  aggregate, 

small  rain. 

2  0311-0318  0°C  large  unmelted  aggregate, 

graupel ,  needles 

3  0318-0320  -  0 . 8  °  C  very  large  aggregate, 

mixed  graupe 1 

4  0325-0327  - 1  . 4  °  C  mixture  aggregate,  needles. 

ice  particles 

5  0327-0338  -2.3°C  more  graupel,  rimed  ice, 

clumps  of  rimed  particles 

6  0340-0359  -4.4°C  needles,  dendrites, 

clumps,  rimed  ice 

o 

0352-descent  unmelted  aggregate  to  4.1  C 


4.6  Summary  of  Case  Studies. 

Figure  4.3  shows  a  summary  of  the  results  of  the  case  studies. 
The  characteristics  of  the  images  are  graphically  depicted  as  a 
function  of  temperature  for  each  of  the  vertical  profiles 
obtained,  and  a  typical  cloud  liquid  water  content  (1  -  30  micron 

droplet  diameter)  is  shown.  This  figure  clearly  shows  a 
systematic  sequence  of  events  as  the  cloud  particles  pass  through 
the  melting  layers.  Those  events  are  as  follows: 

A.  At  temperatures  -2°C  and  colder  therp  is  usually 

a  mixture  of  graupel  embryos  (or  rimed  ice  particles 
of  about  0.3-0. 5  mm),  ice  particles  and  small  aggregates. 
The  graupel  embryos  are  often  associated  with  pockets 
of  liquid  water  content  of  a  few  tenths  of  a  gram  per 
cubic  meter.  These  pockets  are  usually  distinct  from 
pockets  of  ice  particles  and  small  aggregates  where  the 
liquid  water  lias  been  depleted.  The  spatial  dimension  of 
these  pockets  is  on  the  order  of  10  km. 


rv-rvyw t ■  ¥ Tm fvfv." H  x  *Kmjr*rw w. *r.*A*r.>sr  r.- vwm  v* v> 


»T*  .  «  -  ■  .  ' 


KEY 


▲ 

GRAUPEL  EMBRYO 

© 

UMELTED  AGGREGATE 

▲ 

GRAUPEL  PARTICLE 

/ 

NEEDLE 

• 

ICE  PARTICE 

0 

STELLAR 

RIMED  ICE  PARTICLE 

1 

DRIZZLE 

* 

SMALL  AGGREGATE 

• 

VERY  LIGHT  RAIN 

* 

LARGE  AGGREGATE 

•• 

LIGHT  RAIN 

co 

to 

to 

CD 

CO 

0.0  0.2 

03 

CO 

CD 

CO 

1 

CD 

1 

r- 

i 

1 

CO 

1 

1 

CD 

OJ 

1 

r- 

C\J 

1 

N- 

i 

1 

SO 

TYPICAL  CLOUD 

LWC  PROFILE 

CASE 

STUDY 

OATES 

(g  m"3) 

Figure  4.3  Summary  of  data  from  five  case  studies  showing 
vertical  profile  of  cloud  liquid  water  content 
and  evolution  of  hydrometeors. 
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B.  Ice  and  small  aggregates  form  large  aggregates 

beginning  at  about  -2°C.  This  process  proceeds  very 
rapidly  with  large  aggregates  forming  near  the  0°C  level. 

When  present,  pockets  of  graupel  embryos  and  ice  grow 
(slowly  due  to  the  low  LWC)  via  a  wet  accretional  process  to 
1  mm  graupel. 

C.  Near  0°C  there  is  a  mixture  of  large  aggregates 

and  sometimes  graupel  particles,  usually  in  discrete 
pockets  of  ~10  km  across.  Aggregation  appears  to  continue 
through  this  layer  while  the  graupel  particles  start  to  melt 

D.  At  about  0.3°C  small  (1-2  mm)  melting  aggregates  and  large  ! 

(2-5  mm)  unraelted  aggregates  are  observed.  When  present, 
melting  small  (1-2  mm)  graupel  particles  and  1  mm  rain  drops  are 
observed  in  pockets  which  are  usually  separate  from  the 
aggregates . 

o  o 

E.  At  +2  to  +4  C  large  (2-5  mm)  unmelted  aggregates  are 
mixed  with  smaller  (1  mm)  raindrops;  any  graupel 
particles  all  appear  to  have  melted  into  1  mm  rain. 

The  primary  precipitation  mechanism  appears  to  be  through  an 
aggregation  process,  starting  with  0.2  -  0.5  mm  ice  particles 
(often  needles)  near  -3  C  and  small  aggregates,  and  leading 
to  the  formation  of  large  (2-5  mm)  aggregates  near  the  0  C 
level.  These  aggregates  melt  to  form  1  2  mm  rain  with  occasional 
large  (2-5  mm)  unraelted  aggregates  falling  to  "800  m  below  0°C. 

Sometimes,  a  very  different  precipitation  process  is  observed 
which  starts  with  graupel  embryos  ("0.5  ram  rimed  ice  particles)  at 
-3  C,  and  proceeding  through  an  accretional  process  which 
leads  to  the  formation  of  small  (1  mm)  rain  The  aggregation 
process  appears  to  be  more  efficient  in  the  utilization  of  cloud 
liquid  water,  and  we  suspect  that  this  is  due  to  the  long 
residence  time  of  the  (slowly  falling)  aggregates  in  the  region 
above  the  melting  layer,  resulting  in  dry  growth  from  vapor 
deposition.  The  graupel  particles,  on  the  other  hand,  have 
larger  fall  velocities  and  pass  quickly  through  the  region  of 
supercooled  water  above  the  melting  layer,  scavenging  out  a  small 
percentage  of  the  total  volume.  The  relative  importance  of  the 
two  mechanisms  is  the  reverse  of  that  seen  in  summertime 
convective  clouds  where  wet  accretion  is  the  primary 
precipitation  process,  and  aggregation  is  usually  secondary 
(Cooper  and  Lawson,  1984).  This  is  due  to  the  much  higher  initial 
liquid  water  contents  of  convective  clouds  (favoring  accretional 
growth),  and  the  subsequent  rapid  depletion  of  cloud  water  via 
entrainment  (reducing  the  lifetime  shortens  the  residence  time 
for  growth  via  vapor  deposition  and  subsequent  aggregation) 

It  appears  that  the  partially  melted,  large  aggregates,  may 
be  an  explanation  for  the  increase  in  radar  reflectivity.  It  is 
shown  that  these  particles  can  exist  up  to  800m  or  so  below  the 
0  C  level.  The  vertical  profiles  observed  here  are  similar  to 
those  seen  by  other  investigators  (see  Section  1.0  and  Figure 
1.1),  and  it  appears  that  aggregation  may  play  an  important  role 


in  the  formation  of  bright  bands  in  general 


The  curious  part  of  this  chain  of  events  is  how  these  large 
aggregates  develop  so  fast  just  above  the  0  C  level.  The 
"aggregation  zone"  near  0  C  has  often  been  associated  with 
rapid  rates  of  aggregation  (Rogers.  1973).  Rapid  aggregation 
rates  have  also  been  observed  by  Cotton  e ^  a_K  (1986)  at 
much  colder  temperatures  (-10  to  12 0  0  )  in  wintertime  storms 
over  the  Park  Range  in  Colorado.  Numerical  simulations  by  Cotton 
et  al^  (1986)  suggest  that  the  major  term  in  the  aggregation 
process  is  the  scavenging  of  smaller  ice  crystals  by  mature 
(large)  aggregates  via  a  sedimentation  process.  figure  43  shows 
that  at  -3°  to  -4°C  there  are  mostly  single  crystals  and  small 
clumps.  But  at  about  -2°C  aggregation  proceeds  rapidly  through  to 
the  0°C  level  (distance  of  300-400  m)  with  fall  velocities  of 
about  1  m/sec.  That  is.  the  aggregation  process  apparently  took 
only  5-8  minutes,  and  the  rate  of  aggregation  appeared  to 
accelerate  rapidly  near  the  0  C  level,  which  was  also  coincident 
with  the  formation  of  larger  aggregates  which  would  be  required 
to  scavenge  smaller  ice  particles  and  aggregates  via 
sedimentation . 


5.0  EVALUATION  OF  THE  M  METER 


We  have  previously  submitted  a  preliminary  analysis  of  the 
M  meter  in  a  special  technical  report  to  the  A  FOE  Subsequent 
refinement  of  the  data  reduction  procedure  which  includes  a  bear  inc 
friction  factor  and  the  effects  of  heating  of  the  impact  oil 
hydrometeors  was  suggested  by  Mr.  Vern  Plank  of  AFGL,  however, 
lack  of  funds  has  prevented  this  new  analysis  from  being  completed. 
Interested  parties  should  contact  Mr.  Plank  for  a  current 
Interpretation  of  the  performance  of  the  M  meter. 


60  EVALUATION  OF  THE  KA  BAND  RADAR  AND  ATTENUATION 
MEASUREMENTS . 

At  the  beginning  of  the  project  it  was  thought  that  a  nose¬ 
pointing  Ka  band  radar  could  be  used  to  measure  radar 
reflectivity  at  various  ranges  and  from  these  sequential 
measurements  it  would  be  possible  to  deduce  attenuation. 

However,  the  Ka  band  radar  installed  for  this  project  had  a 
higher  noise  level  then  anticipated  and  thus  it  was  difficult  to 
obtain  any  clearcut  measurements  (see  section  3.0  for  examples). 
This  problem  was  amplified  by  the  fact  that  during  the  project  we 
encountered  only  light  to  moderate  precipitation.  The 
reflectivity  returns  were,  therefore,  very  weak  and  often  buried 
in  the  noise. 

The  narrow  beam  and  the  small  pulse  width  led  to  the  sampling 
of  relatively  small  volumes  of  cloud.  From  the  PMS  particle  data 
we  know  there  are  large  variations  in  particle  sizes  over  small 
areas,  hence  significant  reflectivity  variations  would  also  be 
expected  on  this  scale. 

Another  problem  was  that  the  aircraft  altered  its  heading  on 
frequent  occasions.  This  meant  that  the  narrow  beam  was  sampling 
different  cloud  volumes  compared  to  those  subsequently 
penetrated.  Thus  from  a  practical  viewpoint,  only  range  gates 
within  a  few  kilometers  of  the  radar  could  be  used  to  correlate  to 
measured  hydrometeors 


PERSONNEL 


The  following  personnel  worked  on  this  contract: 

C  I  C  Personnel 

Dr.  Larry  G.  Davis  -  Contract  Manager 

R.  Paul  Lawson  -  Principle  Investigator  and  On  Board 

Scientist 

Data  reduction  and  software 
Norman  Ostrander  -  Pilot 
Brian  Jesse  -  Airborne  Technician 

Consultants 

Dennis  Treddenick  -  Detrek  Engineering-Ka  band  radar 

Dr.  Loren  Nelson  -  Ophir  Corporation-software 

Dave  Melita  -  Ophir  Corporation-On  -  board  Scientist; 
reduct  ion 


VMUVVVVWVIV 


8.0  REFERENCES 

Austin,  P  .  M  .  and  A  .  Bern)  s  ,  1050:  A  ijuant  itdtive  study  of  the 

bright  band  in  radar  precipitation  echoes.  J__  Met.eor_x  7, 
165-171  . 

Battan.  L.J.  ,  1  973:  Radar  observations  of  the  a  t_!S  Q  s  g  h  ere_  Univ  . 

of  Chicago  Press,  Chicago.  II.  331  p 

Cohen,  ID.  and  H.J.  Sweeney.  1983  Melting  layer  survey  final 
report.  AFG1.  TR  83  0200,  No.  8-18  .  A! 1 A  H •  1  i 

Cooper,  W .  A . ,  1978:  Cloud  physics  investigations  by  the 

University  of  Wyoming  in  HI  PI,  EX  197  7  Univ  of  Wyoming 
Dept,  of  Atmospheric  Science  R  e  p  u i  f  No  AS  119 

Cooper,  W.  A.  and  R.  P  Lawson,  1981  Physical  interpretation  of 
results  from  the  HIPLEX  1  experiment  3_  Climate  and 
Applied  Me t e o r_x  23,  523  540 

Cotton,  W  .  R  ,  G  .  J  .  Tripoli.  R  .  M  .  Ran  her  a  n  ci  E  A  M  u  1  v  i  h  i  1  1  .  1  986: 

Numerical  simulation  of  the  effects  of  varying  ice  crystal 
nucleation  rates  and  aggregation  processes  on  orographic 
snowfall.  Jx  C_Hma_te  and  Appiped  Meteor_x  25  1658  1  680 

Cunningham,  R.M.,  1947:  A  different  explanation  of  the  bright 

line".  J  x  Meteor^  4.  163. 

Ekpenyong,  B.  and  R.  Srivastava,  1 9  7  o  Radar  characteristics  of 
the  melting  layer,  a  theoretical  study  Univ.  of  Chicago 
Dept,  of  Geophys.  Sciences  and  Illinois  Inst,  of  Tech.  Dept, 
of  Elec.  Eng.,  Technical  Report  No.  16,  Lab  for  Atmos 
Probing,  Chicago,  Illinois. 

Gordon,  G .  and  J . D .  Harwi tz  ,  1  984  An  airborne  comparison  of 

three  PMS  probes.  Jx  Ocean^  l£nbho^xx  _1  ,  2227. 

Hallett,  J.  and  S.C.  Mossop,  1974.  Production  of  secondary  ice 
particles  during  the  riming  process.  Naj_urex  24  9,  2628. 

Knight,  C . ,  1979:  Observations  of  the  morphology  of  melting 

snow.  J _  Atmosx  Sc_ixx  36,  1123  1130. 

K n o 1 1 e n b e r g ,  R . G .  .  1  9  7  6:  The  response  of  optical  array 

spectrometers  to  ice  and  snow:  a  study  of  2D  probe  area  to 
mass  relationships.  Final  Report  Air  Force  Geophysics 
Laboratory,  AFGL  TR  7  6  027  3  M  pp  AD  AO  V- 7  5 1 

Lawson,  R  .  P  .  ,  L  .  G  .  Davis,  and  1)  .  5  T  i  e  d  d  e  ri  i  (  k  ,  1980:  Colorado 

International  Corp.  airborne  seeding  and  cloud  physics 
studies  in  HIPLEX  1977  1980.  Final  Report  to  US.  Dept,  of 

Interior  . 


Lawson,  R.P.  and  R.A.  Stewart.,  1  9  H  :i  An  improved  optical  ice 
particle  counter.  Preprints  Filth  Symposium  on 
Meteorological  Observations  and  Instrumentation,  Toronto, 

5663  pp  . 

Lhermitte,  R.  and  D.  Atlas,  1963:  Doppler  fallspeed  and  particle 

growth  in  stratiform  precipitation.  Proc.  Tenth  Conf.  on 
Radar  Met.,  Washington,  DC.,  297302. 

Marshall,  J.S.  and  K.L.S.  Gunn,  1952:  Measurement  of  snow 

parameters  by  radar.  Jx  Meteor_i  9,  332-327. 

Rogers,  D.C.  ,  1  973:  The  aggregation  of  na  t  uraj_  ice  c  ry  s  t.  a  1_  s_ 

M.S.  thesis  to  the  Univ.  of  Wyoming  Dept,  of  Atmospheric 
Science.  86  pp. 

Ryde,  J.W.,  1946:  The  attenuation  and  radar  echoes  produced  at 

centimetre  wavelengths  by  various  meteorological  phenomena. 
Meteorological  Factors  in  Radio  Wave  Propagation,  London, 

The  Physical  Society,  169-188 

Stewart,  R.E.,  J.D.  Marwitz,  and  J.C.  Pace.  1984:  Characteristics 
through  the  melting  layer  of  stratiform  clouds.  Jx  Atmosx 
Sc  j__x  32,237-273. 

Wexler,  R.,  1955:  Precipitation  growth  in  stratiform  clouds. 

Q  u  a  r  t_  J_  R  oy_  Me  t_e  o  r  x  S  o  c  xx  7  2,  363-37  1. 


APPEND  i  X  A 


PI.  I  GUT  TRACKS 

The  f  o  1  1  o  w  i  n  g  lli  |i  lit  t  i  a  c  k  s  were  u  h  t  a  l  u  e  >1  iis  i  n  it  a 
combination  of  multi  DMF  data,  VOR  I)  MR  data,  and  dead  reckoning 
(integrated  heading  and  true  airspeed)  The  accuracy  of  the 
m  u 1 t i -  DM  E  system  is  about  0.2  km  in  straight  flight  and  0.5  km  l  n 
turns.  On  4-26  and  4  27  86,  the  multi  D M E  w ,i  s  inoperative  and 
the  position  is  plotted  using  only  VOR  DMK  and  dead  reckoning 
The  accuracy  of  the  position  is  degraded  to  about  1  km  in 
straight  flight  and  2-3  km  in  turns.  During  tape  changes  when 
there  are  no  position  data  the  estimated  track  is  shown  in  dashed 
lines.  On  4-7  86,  when  the  aircraft  proceeded  north  of  the 
primary  target  area,  the  position  is  estimated  (but  shown  as  a 
solid  line)  due  to  changes  in  the  cockpit  VOR  and  poor  reception 
from  the  multi-DME.  On  4-26-86  the  aircraft  proceeded  north  of 
the  target  area  in  search  of  precipitataion.  hut  no  usable  data 
were  obtained,  and  thus  the  aircraft  position  from  0232-0330  EST 
is  not  plotted.  The  heading  gyro  did  not  slave  properly  after 
4  8-86  and  the  heading  shown  on  the  48  hr  print  out  (from  the  A  EG  I. 
VAX)  is  in  error.  An  offset  has  been  added  to  correct  these 
headings  on  4-26-86  (135°)  and  4  27  86  (-200°).  The  position 

shown  for  Ha  ns  com  A  F B  is  296.5  true  and  27.6  km,  and  for  the 
Sudbury  radar  is  278.5°  true  and  36.1  km.  It  has  come  to  our 
attention  that  the  postion  that  we  used  for  the  Sudbury  radar  may 
be  slightly  in  error;  if  this  is  the  case,  the  aircraft  tracks 
are  still  valid  and  the  Sudbury  radar  post  i  ion  should  be  adjusted 
on  the  plots. 
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APPENDIX  B 


ITEMIZATION  OF  SOFTWARE  DELIVERED  TO  AFGL 


Program  Name 


£  cri  £_t  _i  on 


AFGL  Reads  9-track  tape,  prints  data  listings 

TAPE  Subroutine  to  read  Bag  tape  fro*  IBM  PC 


AFGL2D 

DRW2D 


Reads  9-track  tape,  draws  PMS  inages 
Subroutine  to  draw  PMS  2d  inages  on  Toshiba 
Printer 


TRKPL 

MLDMT 

MDSR 

SNSR 


Plots  track  from  Multi-DME  position 
Computes  Multi-DME  position 
Subroutine 
Subrout ine 


CPMS  Performs  artifact  rejection  of  PMS  images  and 

computes  liquid  water  content,  equivalent  dBZ 
and  other  statistics  (written  for  HP-1000) 
TWOD  Assembler  subroutine  to  process  image  data 


Performs  spectral  analysis  of  selected 
parameters,  including  de-trending,  tapering, 
power  spectra,  coherence  and  phase  (written 
for  the  HP-1000) 


FFT 


